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^ ■ The chiral transition for two-flavor QCD is predicted to be in the same universality 

^ ■ class as the 3d 0(4) model. This prediction is verified in the Wilson case, but not for 

CN . the staggered- fermion case. The comparison is usually done assuming infinite-volume 

^ i behavior. Here we make an analysis of existing staggered-fermion data using finite-size 

J> ' scaling and normalizing the QCD data. We find better agreement for larger quaik masses. 

^ ■ 

O I In order to check if two systems belong to the same universality class, one looks at 

^ ' universal quantities, such as the critical exponents. More generally, one compares scaling 

T— I ■ functions for both systems. The scaling ansatz implies that the magnetization M of a 

^ ! system be described by a universal function 

^ : M/h'/' = fc{tih'/^') (1) 

where t = {T — Tc)/Tq and h = H/Hq are the reduced temperature and magnetic field, 
^ ■ respectively. Thus, once the non-universal normalization constants To and Hq are deter- 
>■ . mined for a given system in the universality class, its order parameter M scales according 
^ i to the same scaling function fc for all systems in this class. For comparisons with QCD, an 
^ \ important region is the pseudo-critical line, defined by the points where the susceptibility 
X shows a (finite) peak. This corresponds to the rounding of the divergence observed for 
H = 0,T = Tc. The susceptibility scales as x = OM/dH = [1/ Hq) h^l^-^ g{t/h^l^^) 
where g{z) = l/6[fG{z) — z/PfQ^z)]. At each fixed h the peak in x is given by 
tp = Zph^/f^^ and we have Mp = h^^^ fai^p), HqXp = h^^^^^ di^p) ■ Thus, the behavior 
along the pseudo-critical line is determined by the universal constants Zp, fG{zp),g (zp). 

In addition to these infinite-volume scaling laws we may also consider finite-size-scaling 
functions. In fact, the scaling ansatz also implies 

M = L-f^/" Q.ihL^^^/'') (2) 

where L is the linear size of the system and we consider fixed values of the ratio t/h}f^^ = z 
(e.g. z = as in the critical isotherm, or Zp as along the pseudo-critical line). Thus, M 
can be described by a universal finite-size-scaling (FSS) function of one variable. We 
note that in order to recover the infinite-volume expression M = h^^^ fdz) as L — » oo, 
we must have Qz{u) fG{z)u^^^ for large u. Thus, in this limit, the FSS functions 
are given simply in terms of the scaling function fciz)- Working with the FSS functions 
Qz instead of the infinite-volume scaling function fc has the disadvantage that one must 
consider z fixed (thus restricting the regions to be compared in parameter space) but the 
advantage that a comparison can be made already at finite values of L. 
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Figure 1. Comparison at infinite volume. Figure 2. same, using finite-size scaling. 

Critical exponents and the scaling function fc are well-known for the 3d 0(4) model, 
whereas the pseudo-critical line and FSS functions have been recently studied in [1]. We 
find that the peak susceptibilities are given by Zp = 1.33(5) and that the FSS functions 
Qz{u) approach their asymptotic values already at m ~ 2. For QCD, the order parameter 
is given by the chiral condensate < ipip >, the analogue of the magnetic field is the quark 
mass niq, and t is proportional to 6 / g'^ — 6 / g^{0) . Thus, the chiral susceptibility peaks at 
tp ~ m]/^^. Studies at A^^ = 4 by the Bielefeld, JLQCD and MILC groups [2] show that 
the peak positions scale with the predicted exponents, but one sees no agreement with 
the 0{N) scaling function. This comparison was done up to multiplicative constants in 
the fields t and h. We redo it below, performing a normalization of the QCD data. 

Let h = NrrUq/HQ, t = — 6/ g'^{0))/To. From the observed scaling along the 

pseudo-critical line and the universal quantities Zp, fci^p) from the 0(4) model, we deter- 
mine the normalization constants Hq and Tq. This allows an unambiguous comparison of 
the QCD data to fc, as shown in Fig. |I|. We then do a comparison with the FSS functions 
Qz{h L^^^'^) using the data normalized above. We see that along the pseudo-critical line 
(JLQCD and Bielefeld data) FSS works well: the points are already in the asymptotic 
region and scale with the predicted exponent S — as shown in [1] — and with the cor- 
rect multiplicative constant, as expected from our normalization. In fact, all QCD data 
are asymptotic in hL^^^'^, and we can thus rescale all data by their corresponding fciz) 
factors and plot them together, as shown in Fig. ^. We see a better agreement overall, 
although most MILC points are several deviations away from the predicted curve. As in 
Fig. 0, this happens for points outside the pseudo-critical line. Among these points, the 
ones with larger quark masses come closer to the curve. In particular, for the points with 
rUg = 0.025 we see reasonable agreement, suggesting that after finite-size effects have been 
taken into account one may see universal scaling even away from the pseudo-critical line. 
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